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ABSTRACT 

A radio afterglow was detected following the 1998 August 27 giant flare from the soft gamma repeater 
(SGR) 1900+14. Its short-lived behavior is quite different from the radio nebula of SGR 1806-20, but 
very similar to radio afterglows from classic gamma-ray bursts (GRBs). Motivated by this, we attempt 
to explain it with the external shock model as invoked in the standard theory of GRB afterglows. We 
find that the light curve of this radio afterglow is not consistent with the forward shock emission of an 
ultra-relativistic outflow, which is suggested to be responsible for the initial hard spike of the giant flare. 
Nevertheless, shock emission from a mildly or sub-relativistic outflow expanding into the interstellar 
medium could fit the observations. The possible origin for this kind of outflow is discussed, based on 
the magnetar model for SGRs. Furthermore, we suggest that the presence of an ultra-relativistic fireball 
from SGR giant flares could be tested by rapid radio to optical follow-up observations in the future. 

Subject headings: gamma rays: bursts — stars: individual (SGR 1900+14) — ISM: jets and outflows 



1. INTRODUCTION 

Soft gamma repeaters (SGRs) are generally character- 
ized by sporadic and short (~ 0.1s) bursts of hard X-rays 
with luminosities as high as 10^ Eddington luminosity. 
They are also well-know for two giant flares: the first on 
1979 March 5 from SGR 0526-66 (Mazets et al. 1979) and 
the second on 1998 August 27 from SGR 1900+14 (Hurley 
et al. 1999). Frail et al. (1999) reported, following the gi- 
ant August 1998 fiare from SGR 1900+14, the detection of 
a transient radio source. Their observations covered the 
time from about one week to one month after the flare. 
The spectrum between 1 and 10 GHz is well fitted by a 
power-law with F^, cx jy-0-74±o.i5^ rpj^^ source appears to 
have peaked at about a week after the burst and subse- 
quently undergone a power law decay with an exponent of 
a = -2.6+1.5. 

The initial hard spike of the August 27 flare has a dura- 
tion of ~ 0.5s and luminosity greater than 2 x lO^'^^ergs"^ 
(>15KeV) if the source distance is d ~ 7Kpc (Vasisht et 
al. 1994) . The short duration, high luminosity and hard 
spectrum indicate that a relativistically expanding fireball 
was driven from the star. The fireball should be relatively 
clean and the Lorentz factor F ^ 10 was inferred from 
the luminosity and the temporal structure (Thompson & 
Duncan 2001). With the experience of GRB afterglows, 
one may naturally ask whether this power-law fading ra- 
dio afterglow is due to the blast wave emission driven by 
the firebaU. Huang, Dai & Lu (1998) and Eichler (2003) 
had made some discussions on the possible afterglow emis- 
sion from SGRs. In this paper, we try to explain the radio 
afterglow from this fiare. We study the afterglow emis- 
sion from the giant flare in section 2. We find that shock 
emission from an ultra-relativistic outfiow fails to explain 
this radio afterglow. However, we propose that a mildly 
or sub-relativistic outflow expanding into the interstellar 
medium could fit the observations. Finally, we discuss the 
possible origin for this kind of outfiow in section 3. 



2. RADIO AFTERGLOW FROM SGR GIANT FLARES 

We consider that an outfiow with "isotropic" kinetic en- 
ergy Eq and Lorentz factor Fq ejected from the SGR ex- 
pands into the ambient medium with a constant number 
density n. The interaction between the outfiow and the 
surrounding medium is analogous to GRB external shock 
(Rees & Meszaros 1992; Meszaros & Rees 1997), but with 
quite different Eq and Tq. The Sedov time at which the 
shock enters the non-relativistic phase is roughly given by 
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tni — (3i?o/47rnmpC^) — Idays (-Eo,44/"o)^'''j where 
mp is the proton mass and we used the usual notation 
a = 10"a„. As the shock must have entered the non- 
relativistic phase during the observation time of the radio 
afterglow from the giant fiare, we develop a model which 
holds for both the relativistic and non-relativistic phases. 
From the view of the energy conservation, the dynamic 
equation can be approximately simplified as (e.g. Huang, 
Dai & Lu 1999; Wang, Dai & Lu 2003) 
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(F - l)Moc2 + (f2 - l)m,wc' = Eq 



(1) 



where F is the Lorentz factors of the outfiow, mgw = 
{4:/3)TrR^mpn is the mass of the swept-up ISM ( R is the 
shock radius) and Mq is the mass of the original outfiow. 
The first term on the left of the equation is the kinetic 
energy of the outfiow and the second term is the internal 
energy of the shock. 

The kinematic equation of the ejecta is 



dR/dt = /3c/{l-/3), 



(2) 



where v ~ (3c is the bulk velocity of the outfiow with 
^(r) = (1 - F-2)i/2 g^jj^j ^ jg ^jjg observer time. If the out- 
fiow is beamed and sideways expansion with sound speed 
takes place, the expression of mgw and the half opening 
angle of the beamed outfiow 9 are respectively given by 
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where Cg is the sound speed and we use the approximate 
expression derived by Huang, Dai & Lu (2000), which 
holds for both the ultra-relativistic and non-relativistic 
limits. 

Assuming that the distribution of the shock-accelerated 
electrons takes a power-law form with the number density 
given by n{'ye)dje = K-f-Pd-je for 7™ < je < 7m, the vol- 
ume emissivity at the frequency v' in the comoving frame 
of the shocked gas is (Rybicki & Lightman 1979) 
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B^^ KF^ii/y^M, (4) 



where q and rrie are respectively the charge and mass of 
the electron, B±_ is the strength of the component of mag- 
netic field perpendicular to the electron velocity, v'^ and 
v'f^ are the characteristic frequencies for electrons with 
and 7m respectively, and 



F{x)x^P-^^'^dx (5) 



with F{x) = X Kr^i^{t)dt {Kr^j^it) is the Bessel func- 
tion). The physical quantities in the pre-shock and post- 
shock ISM are connected by the jump conditions (Bland- 



ford & Mckee 1976) :n' 
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Y-(r — l)nmpC 



where e' and n' are the energy and the number densities 

of the shocked gas in its comoving frame and 7 is the adia- 
batic index, a simple interpolation of which between ultra- 
relativistic and non-rclativistic limits, 7 = (4r-|-l)/3r, 
gives a valid approximation for trans-relativistic shocks. 

Assuming that shocked electrons and the magnetic field 
acquire constant fractions (eg and e^) of the total shock 
energy, we get 7„ = ee|5f 2i£(r - 1), = VSTreie' 
and K = {p ~ l)n'jP^^ for p > 2. From the spectrum 
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-0.74±0.14 



of the radio afterglow. 



infer that 



F, 

p 2:^ 2.5. It is reasonable to believe that v'j^, in compar- 
ison with the radio frequencies, is very large throughout 
the observations. The observer frequency relates to the 
frequency j^' in the comoving frame by = Di^', where 
D = l/r(l — f3cos6) is the Doppler factor. The observed 
flux density at u is given by 



F, = V^sD^j,,/4nd'' 



(6) 



where Vdi is the effective volume of the post-shock ISM 
from which the radiation is received by the observer and 
should be y = msw/n'nipT^ for the isotropic case. 

Below we shall study the afterglows from the relativis- 
tic outflow with To ~ 10 and the mildly-relativistic one 
respectively. 

2.1. Ultra-relativistic outflow 

In terms of Eqs(l)-(3), we can obtain the dynamic evo- 
lution of the outflow, i.e. we get T{t), R{t), rriswit) and 
9{t). Then using Eq.(4) and the expressions for K 
and 7to, we can get the evolution of the observed flux with 



time. The radio afterglow of August 1998 flare peaks at 
about one week after the burst. In the relativistic shock 
model, it requires that the peak frequency crosses the 
observation band at the peak time if the synchrotron ra- 
diation is optically thin ^. However, this can hardly be 
satisfied for a ultra-relativistic outflow with Fq ~ 10, for 
reasonable values of the shock parameters and the medium 
density n. Instead, the model light curves generally peak 
aX t < O.ldays. This can be clearly seen in Fig. 1, in 
which we plot the model light curves with different values 
of n and ee for both the isotropic and beamed outflow 
cases and compare them with the observation data. More- 
over, the peak flux F^^ are generally much larger than 
the observed peak flux. The reason can be easily under- 
stood from the following analytic estimate for the isotropic 
outflow case. 

At the peak time t ^ lOdays of the radio afterglow, the 
shock had entered the non-relativistic phase, so the radius 
of the shock is roughly R ~ {5/2)f3ct. From the condi- 
tion of energy conservation ^' 



get (3 = {l2E/125TrcH^nmpY^^ 
The magnetic field is -B = 



^-!TR^^nm„c^ = E, one can 
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the peak flux are, respectively,given by 
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where Ne is the total number of the swept-up electrons and 
Pi^,m is the peak spectral power (Sari, Piran & Narayan 
1998). It is clearly seen that Vm can hardly be as large 
as Uohs = 8.46GHz for reasonable shock parameters of Ce 
and eB (e.g. Granot, Piran & Sari 1999; Wijers & Galama 
1999; Panaitescu & Kumar 2002) and, furthermore, the 
peak flux is much larger than detected from the giant flare. 
Though this analytic estimate is for an isotropic outflow 
case, the beamed outflow has also this problem as shown 
in Fig. 1. 

Although the ultra-relativistic shock associated with the 
initial hard spike of the giant flare could not be responsi- 
ble for the observed radio afterglow, we know from Fig. 1 
that its radio afterglow emission should be easily detected 
at the early time even for the beamed case. The optical 
afterglow emission from the ultra-relativistic shock is also 
calculated and shown in Fig. 2. Clearly, early optical after- 
glow emission can be as bright as 100/iJy (R band magni- 
tude niR = 19) at t <; O.ldays for ^ 10~^, a reasonable 
value we know from GRB afterglows. So we urge early 
follow-up radio-to-optical observations for future SGR gi- 
ant flares to test the presence of ultra-relativistic out- 
flows. Even at late time t ^ 10 days, the flux at 
low radio frequencies are intense enough to be detected. 
At, say,;y = 150MHz, the extrapolated flux is 0.4Jy from 
Eqs(7) and (8) for the typical parameters used.^ The X- 
ray afterglow emission from the ultra-relativistic shock is, 
however, predicted to be lower than the detected bright 



^Numerical calculation with the formula about the synchrotron self-absorption process in Wang, Dai & Lu (2000) shows that the synchrotron 
self-absorption optical depth at the peaJj time is several orders of magnitude lower than unity. So, crossing of the self-absorption frequency 
through the observation band is not a viable explanation for this pealc. 

^We estimate that the synchrotron self-absorption frequency is below lO^Hz at this time. 
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and pulsed X-ray afterglow flux (Feroci et al. 2001), which 
is attributed to the emission from the neutron star surface 
immediately after the giant flare. 

2.2. Mildly or sub-relativistic outflow 

In §2.1, we showed that the radio afterglow from the 
1998 August giant flare of SGR1900+14 is inconsistent 

with the shock emission from an ultra-relativistic outflow. 
But a mildly or sub-relativistic forward shock is able to 
explain the observations, as we show below. The reason 
is that for a mildly-relativistic or sub-relativistic outflow, 
it has an enough long period of coasting phase of the out- 
flow, during which the flux increases with time even though 
t'obs- After this coasting phase, the shock deceler- 
ates and the flux begins to fade in a power law manner. 
For an isotropic outflow, the coasting phase lasts about 

(9) 

until the mass of the ISM swept up by the blast wave is 
comparable to the mass of the outflow, where /3o is the 
initial velocity of the outflow. During the coast phase, 
(3 = const and R — /3t oc t^. According to Eqs.(4) and 
(6), (X i?3^(5p-3)/2 ^ ^^^^ ^jjg mildly or sub- 
relativistic outflow is decelerated by the swept-up mass, it 
quickly enters the Sedov phase, during which /? oc t~^^^ 
and R oc t^/s. So, oc t(2i-i5p)/io ^ ^-1.65 ^j. ^ ^ 2.5. 

For a beamed outflow, the coasting phase is similar to 
the isotropic case, as the expansion is dominated by the 
cold ejecta during this phase. But during the deceleration 
phase, the shocked ISM plasma has comparable energy to 
initial energy Eq and the sideways expansion may take 
place, so it is expected that the flux may decay steeper 
than the isotropic case (e.g. Rhoads 1999; Sari, Piran & 
Halpern 1999). 

The fits with model light curves for mildly or sub- 
relativistic outflow are presented in Figs. 3 and 4 for dif- 
ferent shock parameters. In Fig. 3, the isotropic energy 
is chosen to be E^o.iso = lO^^'erg while in Fig 4 the real 
energy of the ejecta is chosen to be lO^^erg (its isotropic 
energy is therefore 10''^/(^|/2)erg, where^j is the initial 
half opening angle). The beamed outflow model can pro- 
vide nice fits of the observations for a wide range of shock 
parameters such as n and Eq. In Fig. 3, we also present the 
model light curve for an isotropic outflow (without side- 
ways expansion) denoted with the dashed line. Clearly, 
it decays too slowly to flt the observations. In all these 
fits, we used fixed values for E, p, 6j, Ce and n with only 
two free parameters: the initial Lorentz factor Tq and e^. 
We therefore conclude that the mildly or sub-relativistic 
outflow from the SGR giant flare could provide a plausible 
explanation for this radio afterglow. 

3. DISCUSSIONS AND CONCLUSIONS 

We have shown that a mildly or sub-relativistic outflow 
from the SGR could be consistent with this radio after- 
glow. This outflow is expected to originate from the neu- 
tron star curst, accompanying the giant flare. SGRs are 
now believed to be "magnetars", neutron stars with sur- 
face fleld of order lO^** — 10^^ Gauss or more (Duncan & 
Thompson 1992; Thompson & Duncan 1995). The en- 



ergy that drives the very large (^ lO^^erg) giant flares 
such as the March 5, 1979 event from SGR 0526-66 and 
the August 27, 1998 event from SGR 1900-M4 is at- 
tributed to a sudden large-scale rearrangement of the mag- 
netic field which releases the magnetic energy, while the 
smaller short repeating bursts {E <^ lO^'^erg) seem to be 
well explained as being driven by the localized yield of 
the crust to the magnetic strain. A magnetic field with 
B > (47r0,„axM)'/' ~ 2 X 10i4(</.„,ax/10-3)i/2G can frac- 
ture the curst, where fi ^ 10'^^(p//9nuc)*^'^crgcm^'^ is the 
shear modulus of the curst and pnuc is the nuclear density 
and (^max is the yield strain of the crust. However, such 
a patch of crust is too heavy to be able to overcome the 
binding energy of the neutron star. We expect that only a 
tiny fraction of the fracturing curst matter can overcome 
the gravitational binding energy and is able to be accel- 
erated to a mildly-relativistic velocity Fq — 1 ~ 0.1 by 
the released magnetic field energy. Note that the kinetic 
energy of the mildly-relativistic matter per unit of mass 
(Fo — l)c^ is comparable to the binding energy GM^s/R, 
where Mns and R are the mass and radius of the neu- 
tron star respectively. Let's denote the amount of matter 
as Am, the isotropic kinetic energy Eq and the real en- 
ergy of the beamed outflow E^ — Eo6j/2, where 6j is 
the beaming angle of this outflow. For Fq — 1 ^ 0.1, 
Am = 5 X 1022E,_43 g = 5 X 10'^^ £0,440"^ g. Let the size 
of this patch of matter be Ar. Because of the insensitivity 
of Ar on Am for the outermost curst of neutron star, we 
estimate Ar ~ 0.1 - 0.3Km for E^ ~ 10-*^ - lO'^'erg. 

Once we know Ar, we can estimate the beaming an- 
gle 6j of the outflow when it breaks away from the con- 
finement of the magnetic field. This amount of mat- 
ter will be vaporized and become plasma near the neu- 
tron star surface, which moves out along the open field 
lines of the magnetar. The initial kinetic energy den- 
sity of this outflow is Sko ~ Er/(A-o/3oc) = 1.6 x 
1024ergcm-3£;^_43(Ar/2Km)-2(/3o/0.4)-i, where E, is the 
real kinetic energy luminosity of this outflow, Pqc is the 
initial velocity of this outflow and Aq is the initial sec- 
tional area. As the plasma moves out to radial radius 
r, the sectional A = 7r(rsin(?)2, where 9 is the an- 
gle relative to the magnetic axis. Because the magnetic 
fleld lines satisfy r oc sin^^, ^ oc r^ for small 9 and 
£fc oc r~^. On the other hand, the magnetic field en- 
ergy density scales with r as = (i?Q/87r)(r/i?)~^, 
where Bq is the surface magnetic field of the neutron 
star. When the magnetic field energy density decreases to 
be comparable to the kinetic energy density, the outflow 
plasma breaks from the confinement of the magnetic field. 
This corresponds to a radius r^/i? ~ (_BQ/87refco)^^'^ — 
30B^/'54;4V'(/3o/0.4)i/3(Ar/0.2Km)2/3 and 9^/00 ^ 

(r/i?)i/2 ~ 5.5Bl%E;ll^{Po/OAy/^iAr/0.2Kmy/^ 
Note that 9i,, which is rouglily equal to the beaming angle 
9j, is very insensitive to the value of Ej-. As the initial 
opening angle near the neutron star surface ^0 = At/R = 
0.02(Ar/0.2Km), we estimate the beaming angle of the 
outflow is 9j ~ 0.1 — 0.2 for typical parameters. 

What powers the ejection of this patch of matter? We 
think that the reconnection of the magnetic field within a 
region of size Ar during the period of the giant flare will 
release energy of {B"^ /8n){ArfVAAt ~ {B^/8n){ArfR, 
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which should be equal to Am(ro — l)c^, where B is the 
crust magnetic field, Va is the internal Alfven velocity and 
At is the growth time of the instability causing the giant 
flare, viz. the duration {At ~ 0.5s) of the initial hard spike 
of the August 27 giant flare(Thompson & Duncan 1995, 
2001). The size is therefore estimated to be Ar ~ 0.2Km 
and the mass is Am ^ g, which are in reasonable 

agreement with the above estimates according to the light 
curve fits if Eq ~ lO^^erg or equally ~ lO^^erg. 

Note that continuing acceleration of electrons at the sur- 
face of the neutron star is also a possible mechanism for 
the radio afterglow, and need further careful investigation 
in future. 

In summary, we studied the afterglow emission from 
the possible ultra-relativistic outflow and mildly or sub- 
relativistic outflow accompanying the SGR giant flares. 
The radio afterglow emission from the August 27 giant 



flare of SGR 1900+14 is consistent with a mildly or sub- 
relativistic outflow, but could not be produced by the 
forward shock emission from an ultra-relativistic outflow. 
However, we predict that this ultra-relativistic outflow, 
suggested to be associated with the hard spike of the gi- 
ant flare, if exist, should produce bright radio to optical 
afterglows at the early phase {t ^ 0.1 days), which can be 
tested by future observations. 
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Fig. 1. — Comparison between the model light curves of the afterglows from ultra-relativistic outflows (Fq = 10) with the observations 
of the radio flare from the August 27 giant flare of SGR 1900+14 observed at the frequency 8.46GHz. Detections and upper limits for the 
non-detections, taken from Frail et al. (1999), are indicated by the filled squares and arrows respectively. The thin solid and dashed lines 
represent the afterglows of isotropic outflows expanding into ISM with n = lcm~^ and n = O.Olcm"^ respectively. Other parameters used are 
^iso = lO^^erg, Ce = 0.3, es = 10~^. The thick solid and dashed lines represent the afterglows of beamed outflows with Oj = 0.15 expanding 
into ISM with n = lcm~^ and n = O.Olcm"^ respectively. Other parameters are the same as thin lines. 
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Fig. 2. — Predicted i?-band {uji = 4.1 X 10^''Hz)optical afterglow light curves of beamed {9j = 0.15) and ultra-relativistic (Fq = 10) outflows 
from SGR giant flares with Ei^a = 10** erg, n = lcm~^ and ee = 0.3, but with different values for Eg. The solid and dashed lines correspond 
to es = 10~^ and cb = 10~^ respectively. 
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Fig. 3. — Fits of the radio (8.46GHz) flare from August 27 SGR giant flare with afterglow emission from beamed, mildly or sub-relativistic 
outflow with isotropic energy Ejso = 10**erg and Oj = 0.15. The solid and dotted lines correspond to beamed outflows with sideways 
expansion expanding into ISM with n = Icm^'^ and n = O.Olcm^'* respectively. Other parameters used are, respectively, (Fq = 1.033, 
Ee = 0.3, es = 0.008) and (Fq = 1.13, ee = 0.3, eg = 0.03). The dashed line have same parameters as the solid line except that no sideways 
expansion is considered. 
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Fig. 4.— Same as Fig. 3 but with Ei^oSy^ = lO^'^'crg. Parameters used arc (n = Icm-^, Fq = 1.13, €e = 0.3, es = 3 X 10"^) for the soUd 
hne and in = O.Olcm"'^, Fq = 1.4, te = 0.3, es = 1.5 x 10~^) for the dotted Une, respectively. 



